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ABSTRACT 

Current semi-analytic models (SAMs) of galaxy formation over-predict the fraction of passive small 
late-type satellite galaxies in dense environments by a factor of two to three. We hypothesize that this 
is due to inaccurate prescriptions on cold gas evolution. In the hope of solving this problem we apply 
detailed prescriptions on the evolution of diffuse hot gases in satellites and on stellar mass loss, both of 
which are critical to model cold gas evolution. We replace the conventional shock-heating motivated 
instant stripping with a realistic gradual prescription based on ram pressure and tidal stripping. We 
also carefully consider stellar mass loss in our model. When both mechanisms are included, the fraction 
of passive late types matches the data much more closely. The satellite over-quenching problem is 
still present in small galaxies in massive haloes, however. In terms of the detectable residual star 
formation rates, gradual diffuse gas stripping appears to be much more important than stellar mass 
loss in our model. The implications of these results and other possibilities, such as redshift-dependent 
merging geometry and tidal disruption, are also discussed. 



1. INTRODUCTION 

In the LCDM model, dark matter structures grow via 
merging and accretion. Within these dark matter halos, 
baryons initially agglomerate onto the centre of the po- 
tential well_witbjii_tbefre^ the h a lo and form 
stars ()Rees fc Ostrikerl HWH IBinnevi HWH ISiil Il977h . 
When the halo becomes large enough, the accreted gas 
can be gravitationally shock-heated and turn into diffuse 
hot halo gas. In the meantime, the constant attraction of 
gravity causes haloes and their galaxies to interact and 
merge, forming clusters and groups of galaxies. While 
orbiting within a host halo, satellite galaxies are likely 
to interact with o ther g alaxies and the hot ambient gas 
(e.g. lChung et^jlll207m lYagi et al.ll2010[ ). The evolution 
of the cold and hot gas contents is heavily influenced by 
the details of these interactions, and the exact mecha- 
nism is not well understood. 

It is known that galaxies in denser environments are 
optically redd er and more quiescent than field galaxies 
(to c it e a few Gisler 119781 : lLarson. Tinslev. fc Caldwell 
Il980t iDresslerl 119801 ) . This strongly suggests that 
the cold gas supply (the source of star formation) 
is controlled by mechanisms closely a ssociated with 
the environment. iGunn fe Gottl (H972) demonstrated 
that the cold gas of a galaxy can be stripped 
off due to the ram pressure exerted on it by 
the intracluster medium as a galaxy moves within 



a cluster potential (see also lAbadi. Moore. fe Bower! 
19991 iQuilis. Moore, fe Bowed 120001: iChung et al.l 120071; 
Tonnesen fe Brvanl 120091 : lYagi et all 120101 ) Tidal strip 



ping is another important process that indirectly affects 
the evolution of cold gas. Satellites trapped by a larger 
cluster halo are tidally stripped off their dark haloes and 
hot gas reservoir. 

Throughout their orbital motions, satellites lose their 
interstellar/diffuse halo gas through shock heating, ram 
pressure and tidal stripping of their dark halo. Conse- 
quently, star formation in satellites may decreases with 



time. Semi-analytic models (SAMs) of galaxy formation 
assume that the shock heating is very efficient, and it is 
generally assumed to ins tantly remove diffuse h alo gas 
from satellite haloes (e.g. Tva uffmann et al1l!999h . Since 
the hot gas reservoir is removed instantly, gas cooling 
stops, cold gas is quickly depleted through star forma- 
tion and SN feedback, and eventually star formation is 
quenched. Such models predict that the bulk of satellite 
galaxies in large haloes should be virtually passive, with 
only about 20% predicted to be active (e.g. iKimm et al.l 
2009, hereafter K09). 

However, this prediction is not supported by ob- 
servation. Studies based on the specific star for- 
mation rate s measured f rom GAL EX and SDSS ob- 
servations (jSalim et al.1 12007! ) have found that a 
much larger frac t ion o f satellite galaxies are ac- 
tive (jKimm et al.l l2009f) . This is sometimes re- 
ferred to as the sate l lite over- quenching problem (e.g. 



Weinmann et al.l20 06; Baldrv et al. 2006; van den Bosch 
20081: iGilbank fe Baloghl 120081 : iFontanot et al.l 120091 : 
Kimm et alll2009l ). 



The assumption of instantaneous shock heating of 
satellite gases is thought to be the most likely cause of 
this mismatch between model and data. Observationally, 
there is evidence for hot gas around satellite galaxies, 
as shown by Chandra observations of ea r ly- type galax- 
ies in nearby clusters (jSun et aT~l 120071 : Ueltema et al.1 
2007). These observations are bolstered by recent hydro- 
dynamical simulations suggesting that a non-negligible 
fraction (up to ~ 40%) of diffuse hot gas remains associ- 
ated with satellite galaxi es for several Gyrs after merg- 
ing with the host halo ([Bekki. Couch, fe Shioval 120021 : 



iKawata fc Mulchaevll2008HMcCarthv et al.ll2008D . Vari- 
ous S AMs have included the effect of ram pressure strip- 
ping (iKhochfar fc Ostri kcr 2008; Ka ng fc van den BoschI 
l2008t IFont et al.l 120081: IGuo et al.ll2010l ). In a recent at- 
tempt iFont" et al. n^ooar showed that inclusion of ram 
pressure stripping can help improve the bimodal colour 
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distribution of satellite galaxies and the dependence 
of satellite colours on gala xy environment (see also 
iKang fc van de n Bosch 120081 ). However, they did not 
take into account the decay of the satellite orbits due 
to dynamical friction, and hence their calculations may 
slightly overestimate the amount of hot gas retained in 
the satellite systems. Further, as we will show, the in- 
stantaneous recycling approximation for stellar mass loss 
used in the model may change the r ecent star forma- 
tion h istory of galaxies. Alternatively, iWeinmann et al.l 
propose a simple gradual recipe for the stripping 
of gas and dark matter that reproduces the specific star 
formation rates as a function of clustocentric radius. 

Another factor that can contrib ute to the interstel- 
lar m edium is stellar mass loss (e.g. iBregman fc Parriotti 
2009) . Stellar masses vary widely and their lifetimes de- 
pend strongly on their mass. Thus, the stellar mass loss 
of a population is not instantaneous but changes grad- 
ually with time (though the overall change can be dra- 
matic). Yet, most SAMs (|Hatton et a l. 2003 is a notable 
exception) have used a simple approach in which roughly 
30% - 40% of newly formed stars insta ntaneously evolve 
off a nd are recycled to cold gas (e.g. iSomerville et al.l 
2008) . This may be a reasonable approximation for long- 
timescale star formation, but it may not be adequate for 
investigating short-timescale phenomena, such as the re- 
cent star formation history of satellite galaxies. 

In this study, we aim to investigate the effects of the 
two physical processes mentioned above in the context of 
SAMs. We first consider the orbital motion of a satellite 
under the influence of dynamical friction and the associ- 
ated tidal stripping, both of which have been neglected 
in th e previous studies based on EPS (e.g lLacev fc Cold 
f 993) formalism^. We also include a prescription for ram 
pressure stripping (§2.1). Next, in §2.2, we present fur- 
ther improvements on the earlier implementations of stel- 
lar mass loss. We finally describe the impact of these con- 
siderations on satellite galaxy evolution in §3 and discuss 
the implications of this in §4. 

2. SEMI- ANALYTIC MODEL 

We adopt a semi-analytic approach of galaxy formation 
to investigate the effect of different physical processes 
on the recent star formation history of satellite galaxies. 
The fiducial SAM w e are u s ing is described in deta il in 
IKhochfar fc Burkerti (120051) : IKhochfar fc Silk) (f2006h . In 
what follows we briefly lay out the basic ingredients of 
the SAM and explain the important modifications to the 
fiducial model. 

We generate dark matter halo merger trees for a cosmo- 
logical volum e of 10 6 Mpc 3 using a Mont e- Carlo method 
developed by ISomerville fc Kolattl (|1999[ ). Within these 
halos gas cools and collapses at the centre of each halo, 
forming a rotationally-supported disc galaxy. Star for- 
mation in this disc is modeled assuming a Schmidt- 
Kennicutt law (|Kennicuttl 119981 ) . Some massive stars 
quickly die as supernovae and release their energy, which 
heats up the cold gas in the disk. Note that our fidu- 
cial model does not include recycling of mass loss by 

1 SAMs based on N-body merger trees can follow the tidal strip- 
ping of dark matter subhaloes as long as they are robustly iden- 
tified. Beyond that point, satellites are tracked based on analytic 
formalism. 



stellar evolution. When two haloes merge, the most 
massive galaxy in the more massive halo is set to be 
the "central" while the other galaxies become "satel- 
lite galaxies". These satellites begin to fall toward the 
central galaxy due to dynamical friction and eventu- 
ally merge with the central galaxy. If the mass ratio 
(= Mhost/Msat) between two galaxies is less than 3.5, 
we assume that the existing stars in both galaxies and 
the new stars formed by merger-driven starbursts build 
up the spheroidal component of the merger remnant. It 
should be noted that galaxy m ergers are the only w ay to 
form a bulge in our model (c.f lAthanassou lal l2008l ). We 
use a mass ratio-dependent burst efficiency to determine 
the number of stars formed duri ng the merger, d rawn 
from hydrodynamic simulations (|Cox et al.l [2008) . As 
galaxies grow, black holes develop as a result of galaxy 
mergers (Kauffma nn fc HaehneltlT2000[ ) and become ac- 



tive in suppressing star formation in their host galaxies 
(jSchawinski et al.1 12006'). Throughout this paper, we use 
the following set of cosmological parameters: fio = 0.27, 
A = 0.73, n h /n = 0.15, ct 8 = 0.77, and h = 0.71. 

2.1. Dynamical Friction and Gas Stripping 

Earlier implementations of the dy namical friction 
of sat ellite haloes were based on the [Chandra sekharl 
(1 19431) formula with a fixed satellite halo mass (e.g. 



Kauffman n et al.l 119991 see however iTavlor fc~ Babul 
2004T for models with mass loss). Detailed nu- 
merical simulations, however, have found that 
the Chandrasekhar formula systematically un- 
derestimates the merging timescales, especially 
for large mass ratios (Mh os t/M sat ^> 1) (e.g. 
iBovlan-Kolchin. Ma. fc Quataertl(2008l) . 

To delineate the evolution of satellite systems more 
precisely during halo mergers, we calculate the two- 
dimensional orbital motions and tidal mass-loss of the 
satellite halo. We treat each subhalo as a point-like ob- 
ject orbiting in an isothermal host potential. Satellite 
galaxies embedded in a subhalo are initially placed at the 
virial radius of the host halo at the beginning of merg- 
ers, orbiting with the circular velocity of the host halo 
(jBensonl 120051 : IKhochfar fc Burkerti I2006Q . Tangential 
and radial velocities of the infalling haloes are assigned 
as (r]V c , y/l — Tj 2 V c ), where n is the circularity defined as 
the ratio of the orbit angular momentum to that of cir- 
cular orbit with the same energy. We ado pt the random 
circularity following ILacev fc Cold ( 19931). The exerted 
dyna mical friction is computed as (Binncv fc Tremaind 
1998). 



dv 
ctf 



GM SSLt {t) 



In A 



(1) 



erf 



exp 



(2) 

where A is a Coulomb logarit hm and we use the v alue 
of (1 + A/host /M 8at ) following ISpTingel et alJ (120011 ). V c 
is the circular velocity of the halo at the virial radius, v 
is the orbital velocity of the satellite halo, e is the unit 
velocity vector and r is the distance of the satellite from 
the centre of the host halo. Here the satellite mass M sat 
explicitly depends on time. 
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To evaluate the satellite mass M sat during its orbit, we 
calculate the radius within which the motion of a particle 
is governed by the satellite. We assume that all satellite 
mass within the sphere of influence (r so ; ) is bound to the 
satellite, while the matter outside of r so j is strippe d. The 
sphere -of- influence radius (r so i) can be written as (|Battinl 
1987) 



Msat (t) 

M host {< r) 



-0.4 



(1 + 3 cos 2 



0.4 cos 9 



\0.1 



1 + 6 cos 2 9 
1 + 3 cos 2 e 



(3) 



where r is the separation between the central and satel- 
lite galaxies, M sat (i) is the total (baryon+dark matter) 
mass of the satellite halo, Mh os t(< r) is the total mass 
of the central halo within r, and 9 is the angle between 
the line connecting the particle to the center of mass 
of the satellite halo and the line connecting the centers 
of mass of the satellite and the host haloes. For large 
mass ratios (Mh os t/M sat >• 1), Eqn. [3] can be reduced 
to r soi ~ r[M sat (t) / M hos t(< r)} 2 / 5 . At each time step 
(St), we assume that a fraction (St/td yn ) of dark mat- 
ter outside r SO j is stripped, where ta yn is the dynamical 
time scale of the satellite halo. We assume that a galaxy 
merger takes place if the separation between the two 
haloes becomes smaller than the size of a central galaxy, 
which is computed from the ob servational relati onship 
between galaxy mass and size (IShen et all [20031. We 
assume the size of a galaxy to be twice the effective ra- 
dius. By doing so, the merging time scale obtained from 
our simple analytic m odel shows good agreement with 
the nu merical study of lBovlan-Kolchin. Ma, fc Quataertl 
(2008) (within 20% error ) . 

It is interesting to compare our simple model for mass 
loss with other studies. Our calculation gives a mass loss 
of approximately 40% for the dark matter component 
during the first pe ri-centric passage , whic h is comparable 
to what is found in lTavlor fc Babull (|2004T ) . Note that our 
calculation does not include the effect of tidal heating. 
We have also compared our mass-loss curves with results 
obtained through N-body simulations using GADGET-2 
(Yi et al. in prep), and have found qualitatively good 
agreement. 

As a satellite orbits in the host potential, the diffuse 
gas of the satellite will be stripped due to the ram pres- 
sure exerted by the intracluster medium ()Gunn fc Gottl 
[19721) . We mo del the ram pressure fol lowing t he analytic 
formu lation oflMcCart hv et ahl (|2008[ see also lFont et al.l 
2008) . We first compute a radius r ram beyond which ram 
pressure is able to strip material 



Phost(r)v 2 lh (t) 



'2Cji A 1 (<C ?"ram)Pgas(^ram) 



(4) 



where phost is the density of the host halo gas, w or b is the 
orbiting velocity of a satellite, and p gas is the gas density 
of a satellite halo. Similar to r soi , we remove a fraction 
(dt/t CTOSS ) of diffuse gas outside r ram . However, if r so ; is 
smaller than r ram , we take r so ; as the stripping radius. 
Using this analytic approach, we obtain halo gas loss of 
approximately 60% after the first pericentre passage for 



(M host /M sat , v r /v c , vt/vc) = (25,0.9,0.7), which is com- 
parable to the results of McCarthy et al.l (|2008ft . It is 
worth me ntioning that thi s gas stripping prescription is 
similar to lGuo et al.l (|2010f ) . however the two models dif- 
fer slightly in several manners. For example, they adopt 
a fixed orbiting velocity for satellite haloes, whereas our 
model follows the evolution of the orbital motion until 
the galaxy at the center of the halo merges with the cen- 
tral galaxy. 

As the outer part of the satellite diffuse halo gas is 
stripped during its orbit, we assume that the remain- 
ing gas quickly restores an isothermal profile by redis- 
tributing the mass inside the stripping radius. Thus, 
the cooling rate becomes smaller as more diffuse gas is 
lost. Cold gas replenished by the satellite halo gas can 
subsequently form stars, and supernova explosions are 
assumed to blow cold gas out into the satellite diffuse 
halo. This assumption affects the star formation history 
of satellite galaxies, and we will discuss the implications 
in §4. 

2.2. Stellar Mass Loss 

A large fraction of stellar mass is believed to be re- 
leased in the form of cold gas and becomes available 
as a source of the next generation stars. Stellar evo- 
lution and mass loss strongly depend on stellar mass it- 
self. Low-mass stars evolve with a long characteristic 
timescale, whereas high-mass stars are short-lived and 
end with dramatic mass loss processes. This difference 
gives an important characteristic of stellar mass loss for 
each stellar population. Since the mass spectrum of the 
population is continuous, stars evolve off the population 
either through supernova explosion or a planetary nebula 
process continually. 

Conventional SAM models often neglect this fact be- 
cause tracing stellar masses from the entire stellar sub- 
populations is a computationally expensive task. In- 
stead, they generally calculate the fraction of ejecta for a 
given stellar population and apply it to all other popula- 
tions regardless of age. For exa mple, when the Scalo ini- 
tial mass function (Scalo 1986) is adopted, a stellar pop- 
ulation will return roughly 30 percent of its mass during 
its evolution. Then, the models assume that 30 percent 
of newly formed stars quickly evolve off and constitute 
the galaxy's gas component. This may be a reasonable 
approximation for studies of the long-timescale star for- 
mation history of galaxies, but it is not as applicable to 
short-timescale histories, which are the main focus of our 
investigation. 

In this study we develop an economic prescription 
whereby stellar mass loss can be approximated with 
reasonable accuracy. Since massive young stars evolve 
off quickly and their remnant fraction in mass is low 
(i.e., a higher mass loss rate), a large mass- loss rate 
can be expected for galaxies with young populations. 
On the other hand, if the mean age of a galaxy is 
high, low-mass stars dominate the stellar mass loss, 
thereby lowering it. Therefore, one can anticipate the 
anti-correlation between the stellar mass- weighted age of 
a galaxy ((t) m ) and the specific stellar mass- loss rate 
( dTO [ oss = dmi osa /dt/m sa i). 

In order to obtain a relation between dm[ oss and (t) , 
we fi rst compute the lifetim e of stars with a broken-power 
law ([Ferreras fc Silld l2000f ) . which is obtained from the 
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data of lTinslevI (fl980h and ISchaller etaD (fl99l . 

Z™_ = { 9.694 (m/M )" 2 - 762 m < 1OM , , 
Gyr 1 0.095 (m/A/©)" 764 m > 10M Q . [ °> 

At the end of their lifetimes, stars will release most of 
their mass, leaving small remnants such as white dwarf 
or black holes. The m ass locked up in the r emnant is 
approximated following iFerreras fc Silkl (|200Q[ ) as: 

f O.l(m/M ) + 0.45 to < 10M Q 
^ = \ 1.5 10 < m < 25M 

M © { O.61(m/M ) - 13.75 to > 25M . 

(6) 

Using this, we compute the stellar mass loss of galaxies 
with a variety of star formation histories, drawn from 
the fiducial SAM without our new features. Assuming 
that newly formed stars at each redshift follow the Scalo 
initial mass function, the specific stellar mass-loss rate 
is computed as a function of (t) m following cosmic star 
formation histories. This (i) m — dm[ oss relation will be 
used to estimate the approximate amount of stellar mass 
loss for a given simulation period (= dTO[ oss TO g aidi) in our 
SAM. This approximation is very reliable for ordinary 
star formation episodes, but it underestimates the mass 
loss when star formation rates are high (e.g., if there is 
a merger-driven starburst in the previous time step). In 
order to compensate for this shortcoming, we explicitly 
add the mass loss from the youngest stellar population 
in the model galaxy to the mass loss estimated from the 
(t) — ^ m Ioss relation. This treatment can be expressed 
as follows: 

(5toi oss ~ [dm[ oss (t m ) TO gal + dm' loss (t y ) Sm y ] St (7) 

where t m is the mass-weighted mean stellar age of the 
target galaxy, t y is the mass-weighted mean age of 
its youngest populations, St is the integration timestep 
used in our SAM («250 Myr at z = and «80 
Myr at z = 1), and Sm y is the mass of the stars 
formed since the last time step. Our modeling of non- 
instantaneous stellar mass loss is sim pler than the re- 
cent m odels of lArrigoni et a 1 H2010D or IBenson fc Bowerl 
(2010). They fully consider the mass loss for the past 
time-s teps, and include a metallicity-dependent stellar 
yield (|Portiriari. Chiosi. fc Bressanl 119981: iMarigol I2001D 
in the latter paper, which is neglected in this work. 

Fig. [T] shows the specific mass-loss rates of galaxies 
having a wide variety of star formation histories as grey 
contours and the (t) m — dm[ oas relation as a solid line. 
The specific mass-loss rate of a simple stellar population 
is also shown as a dotted line. Since subpopulations of a 
galaxy have different ages and younger populations yield 
more mass loss, the total dm[ QSS is usually greater than 
that of a simple stellar population. As expected, the 
specific stellar mass-loss rate is a strong function of the 
mass- weighted age of a galaxy. It is interesting to note 
that stellar mass loss rates are so small that only a popu- 
lation younger than lOOMyr can supply enough cold gas 
to make galaxies switch from passive to active, assuming 
a star formation efficiency of 2% (Eq. [5]). 

A major source of uncertainty in our (t) m — dm[ oss 
approximation is the fraction of stellar mass loss that 
contributes to the diffuse halo gas and cold disc gas. 
iParriott fc Bregnianl (|2008l ) show using two dimensional 
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Fig. 1. — The relation between the stellar mass- weighted age 
of galaxies ((t) m ) and specific stellar mass-loss rate (dmJ oss = 
dm\ oss / At / m && \) . Stellar mass-loss is estimated for galaxies with 
a variety of star formation histories, drawn from a semi-analytic 
model (gray shading). The mean and lc error of this relation is 
shown as a solid line. The specific stellar mass-loss rate obtained 
from a simple stellar population with the Scalo IMF is shown as 
a dashed line. The (t) m — ^ m J osa relation is used to approximate 
the continuous stellar mass loss at each time step. 

hydrodynamic simulations including radiative cooling 
that 25% of the mass loss from red giants remains cool. 
On the other hand, most of the mass loss from the 
planetary nebula phase app ears to remain warm or cool 
(jBregman fc Parriottll2 009) . Since the main contribution 
of the mass loss comes from stars with masses of 2 - 8 
Mq, whose final stages are planetary nebula, we assume 
that half of the mass loss goes into the hot halo gas and 
the rest returns to the cold disc gas reservoir. The exact 
values of the fractions are unknown, but the models are 
not very sensitive to such variations (within a factor of 
two). 

2.3. The effects of the new prescriptions 

In Fig. [5] we show the impact of the two different 
physical prescriptions on the cold gas evolution for a 
2 x 10 10 Af Q disc-dominated satellite galaxy. The satel- 
lite galaxy can retain a copious amount of cold gas from 
diffuse gas cooling if it maintains a hot gas halo for an 
extended period of time. This leads to a further growth 
in stellar mass in the satellite galaxy. It should be noted 
that the contribution from stellar mass loss is an order 
of magnitude smaller than the prolonged gas cooling in 
this specific example, and results in an order of magni- 
tude increase in cold gas at z = 0, compared with our 
fiducial model where no recycling is used. 

3. RESULTS 

In this section, we investigate the effect of the two 
aforementioned physical processes on the recent star 
formation history of satellite galaxies. There are var- 
ious ways to quantify recent star formation activ- 
ity; we make use of the specific star formation rate 
(SSFR=SFR/TO g ai), obtained from UV and optical pho- 
tometry. UV lights are much more sensitive to current 
and recent star formation than optical colours, and, com- 
pared with emission line diagnostics, they detect not 
only on-going but also recent (up to about one billion 
years) star formation. Empirical specific star forma- 
tion rates are the values averaged over the past 100 Myr 
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Fig. 2. — The impact of the different physical prescriptions on 
the cold gas evolution in a sample satellite galaxy. The vertical 
dashed line marks the moment where a central galaxy becomes a 
satellite galaxy through a halo merger. The fiducial model (solid 
line) assumes that the hot gas halo of the satellite is shock heated 
and instantly stripped away at the time of the merger. Note that 
the fiducial model does not take into account stellar mass loss. The 
dotted line shows the model in which stellar mass loss contributes 
to the cold gas disc. The dashed line shows the model that includes 
the gradual diffuse gas stripping prescription but not the stellar 
mass loss. The dash-dot line shows our new model, which includes 
both features. The rapid drop in cold gas around 8 Gyr is due to 
a major merger that turned much of the cold gas into stars. 

(jSalim et al.1 120071 ) and model values are averaged over 
the past 270 Myr. 

We consider a galaxy "passive" if the following condi- 
tion is met: 

logSSFR < -9 - 0.21og(m gal //i- 2 Af Q ), (8) 

We choose this threshold because GALEX ultraviolet de- 
tections and star formation rate measurements are robust 
above this level. This criterion reasonably divides galax- 
ies into active and passive galaxies (see K09). 

The observational sample used here is the same used 
in K09, which is drawn f rom the GAL EX-SDSS matched 
sample constructed by iSalim et al.l (l2007[). K 09 also 
make use of the Group catalog |Yang et aLll2007l ) for es- 
timates of the dark matter halo mass and galaxy stellar 
mass. For a more detailed description of the data, read- 
ers are referred to K09. 

3.1. Recent Star Formation And Morphology 

IKimm et al.l (|2009D inspected star formation rates with 
respect to galaxy mass and halo mass, but did not 
consider galaxy morphology. It is important to check 
whether the satellite over-quenching problem is present 
in all morphological types of galaxies. In order to 
examine this, we divide the observational data into 
early and late types using the concentration index in 
the r-band by r equiring that early t ypes have C r = 
jW^ sn > 2.6 (|Strateva et al.l 120011 : iShimasaku et al.l 
2001)). Model galaxies are classified based on stellar 
mass- weighted bulge-to-total (B/T) ratio. We assume 
that early types have a B/T ratio of greater than 0.4 
(Fig. Et, which is motivated by o bservational studies 
(e.g. lSimien fc de Vaucouleurslll986l ). where the division 
occurs around B/T=0.4. 

Fig. [4]shows the fraction of passive satellites (/pass) for 
different morphologies. The observational data (Fig. [4j 
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Fig. 3. — The fraction of early-type galaxies (/ ea rly) as a function 
of galaxy stellar mass. We assume that early-type galaxies have 
r-band concentration indices greater than 2.6 in the observational 
sample. 



(a)) show a strong mass dependence of the passive galaxy 
fraction. Interestingly, the mass dependence of / paS s vir- 
tually disappears when focusing on individual morpho- 
logical classes (Figs. 0J-(b) and (c)), which implies that 
the overall mass dependence is a result of a combination 
of scaling relations among g alaxy mass, type , and star 
formation rates (Fig. g} (e.g. IPark et al.ll2007D . 

The morphological mixture as a function of galaxy 
mass is reasonably reproduced by our fiducial SAM 
(Fig. 0}. The passive galaxy fractions for early types 
are reasonably reproduced by most models as long 
as they include A GN feedback (jSchawinski et al.l I2006I : 
IKimm et~a l. 2009) . All models present here include AGN 
feedback. However, most late-type model galaxies do 
not match observations, as they generally predict that at 
least 80% of late types must be passive (Fig. H}(c)). As 
discussed previously, late-type galaxies cannot maintain 
the observed level of star formation, most likely because 
of the lack of cold gas. 

Fig. EJ-(c) shows how our new prescriptions change the 
model outputs against the observed late-type galaxies. 
Model galaxies become progressively more active when 
stellar mass loss and gradual gas stripping are considered. 
The inclusion of stellar mass loss alone still produces too 
many small passive spirals. This is because the amount 
of stellar mass loss released from old stellar populations 
in small galaxies is too little to fuel star formation. When 
both gradual diffuse gas stripping and continuous stellar 
mass loss are adopted, the passive galaxy fraction be- 
comes closer to the data. As an indication of the sensi- 
tivity of the observational fraction of passive galaxies to 
the actual concentration index cut C r , we choose a cut of 
C' r = 2.4 as an example. This cut is more conservative in 
the sense that it provides a cleaner late-type sample. We 
find that the more strict cut results in a lower passive 
galaxy fraction (/ paS s ~ 0.1 — 0.3), indicating possible 
contamination from early-type galaxies within the sam- 
ple selected using the less conservative cut. It should be 
noted that our models predict a strong dependence of the 
passive fraction on galaxy mass for late type satellites, 
which is not observed. This indicates that the evolu- 
tion of cold gas in late type satellites is not yet correctly 
captured by the model. 
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Fig. 4. — The fraction of passive galaxies (/ pa aa) as a function 
of galaxy stellar mass. The top panel shows the / paS s for all satel- 
lite galaxies. Early-type (middle) and late-type (bottom) galaxies 
are also shown separately. The empirical data is indicated by solid 
lines, while our fiducial model (old), the model with gradual diffuse 
gas stripping (GDS), the model with stellar mass loss (SML), and 
the model with both GDS and SML are shown by black dotted, 
green dashed, orange dottcd-dashed, and red solid lines, respec- 
tively (see chart key). Our original semi-analytic model overpro- 
duces passive late types. Note that late- type galaxies become more 
active (consistent with observations) when GDS and SML are take 
into account. 



The inclusion of stellar mass loss and gradual diffuse 
gas stripping affects the morphology of model galaxies in 
a counterintuitive way. Fig. [3] shows that such consid- 
erations raise the early-type galaxy fraction. This may 
sound counterintuitive, because those processes gener- 
ally enhance the star formation in the disc. However, 
the morphology (the bulge-to-total ratio) of a satellite 
galaxy rarely changes due to these considerations. The 
increase in the disc mass due to prolonged star forma- 
tion can be large enough to make a galaxy "active" , but 
is too small to change the galaxy morphology. On the 
contrary, an increase in the stellar disc occurring in the 
central galaxy leads to increases in the bulge mass of the 
central galaxy whenever it merges with another galaxy, 
according to the current semi-analytic prescription. This 
has the net effect of increasing the number of galaxies 
with a higher value of bulge-to-total ratio. 

Our criterion for passive galaxies is rather conservative, 
and it is interesting to note that relaxing this criterion 
to include intermediate colors alleviates the discrepancy 



between model and observation. This implies our model 
produces too many inte rmediate-color late types (see also 
iWeinmann eTalll20Tol ). 

3.2. Environmental Dependence 

The observed fractions of passive satellite galaxies for 
different environments (halo mass) and morphologies are 
shown in Fig. [5] (top row). The estimates of dark matter 
halo masses for observed galax ies are adopted fro m the 
Group catalog constructed by lYang et all (|2007t ) . The 
colour-codings represent different halo masses. As can 
be seen, both early and late types show an environmen- 
tal dependence in that star formation activity in larger 
clusters is suppressed. Again, when the morphology is 
fixed, the dependence of / paS s on galaxy mass and halo 
mass is weak. 

The environmental dependence can also be found in 
our original (fiducial) SAM (the second row from the 
top) where the diffuse halo gas of satellites is instantly 
shock-heated at the beginning of halo mergers. This can 
be understood by noting that satellite galaxies in more 
massive halos merge with their host halo earlier and use 
up their available cold gas without any replenishment. 
This mechanism is even clearer for late types, because 
early-type satellites are m ore strongly influenced b y AGN 
feedback in our models (Scha winski et al.l |2006( ). Note 
that more than 80% of the small late types in the smallest 
haloes (black line) are predicted to be passive in this 
version, which is inconsistent with observation. 

Models with gradual gas stripping (third row from the 
top) shows a better match to the data. The extended 
cold gas supply preferentially lowers the passive galaxy 
fraction of massive galaxies. This is because the gravi- 
tational restoring force of such systems is usually strong 
enough to retain a large amount of halo gas. Since the 
ram pressure becomes less effective if the mass ratios 
between the host and the satellite halo are comparable, 
galaxies in less massive haloes are more likely to be ac- 
tive, for a given galaxy mass. However, this model still 
predicts that a substantial fraction of small late types in 
massive haloes are passive. Hence, we argue that gradual 
gas stripping alone cannot reproduce the observed level 
of recent star formation activity. 

The effect of stellar mass loss on recent star formation 
history appears to be smaller than that of gradual diffuse 
gas stripping (fourth row from the top). Most galaxies 
in large haloes are predicted to be passive, which is in- 
consistent with observation. It seems that stellar mass 
loss alone cannot maintain the observed level of star for- 
mation of late types in cluster environments. This is still 
true even if we arbitrarily allow 100% of the stellar mass 
loss to contribute to the cold phase gas. Lastly, the model 
with both gradual diffuse gas stripping and stellar mass 
loss shows a further improvement in matching the ob- 
served fraction of passive small satellites (bottom row). 
Since cold gas supplies are established through two chan- 
nels, the overall passive galaxy fraction is diminished. 
This model slightly underproduces the massive and pas- 
sive late types. When the total sample is considered (left 
column), the models appear to reproduce the observed 
dependence of passive galaxy fraction on galaxy stellar 
mass. Based on our simple exercise using gradual dif- 
fuse gas stripping and stellar mass loss, we conclude that 
both contribute to the recent star formation history of 
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Fig. 5. — The passive fraction of satellite galaxies (/ paS s) as a function of galaxy stellar mass, for different halo mass bins, as indicated 
by different colours and symbol sizes (see plot legend). Top panels display the observed values. The theoretical models with different 
prescriptions are shown in the subsequent rows. We present the results for different morphologies, as indicated in the plot. The SAM 
with gradual diffuse gas stripping (GDS) and stellar mass loss (SML) shows improved agreement with the data, although small galaxies in 
massive haloes are admittedly still too passive. 



satellite galaxies. 

4. LIMITATIONS 

We have shown that allowing a diffuse gas reservoir and 
continuous stellar mass loss for satellite galaxies can fur- 
ther supply cold gas, and hence may increase the fraction 
of actively star-forming spirals in group environments. 
However, we note that the amount of the supplied gas 
that is actually used to form stars depends on the mech- 
anism of the supernova feedback. 

As an extreme case, if supernovae do not strongly blow 
out cold gas from discs, more would be available for star 
formation, and the fraction of passive spirals would de- 
crease. Additionally, satellite galaxies would form stars 



for a long period of time (e . g lKennicutd[l998h even if dif- 
fuse gas reservoir were depleted. On the other hand, if 
supernova explosion blow most of the cold gas away from 
the satellite halo, there would be less hot gas remaining. 
Since the cooling rate is comparable to the blowout rate 
by supernovae in small subhaloes, where the disagree- 
ment in the passive fraction is the most notable, the 
amount of cold gas would drop rapidly. Then, most of 
the small galaxies that have already moved around their 
central galaxy several times are more likely to be passive 
at z = even if they may have their own hot gas halo. 
This may suggest that satellite over-quenching might not 
be solely due to inadequate treatment of environmental 
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effects, but may also relate to how one models supernova 
explosions. 

In this study, we have assumed that supernova explo- 
sion cannot blow cold gas away but instead transforms 
it to hot diffuse gas. Clearly, this is an oversimplifica- 
tion, and further investigation is needed on this subject. 
It is interesting to note, however, that relaxing the as- 
sumption so that the reheated gas can escape the satellite 
potential would exacerbate the mismatch in the passive 
fraction for small disc-dominated satellite galaxies. It is 
also possible that the supernova feedback is ineffective in 
transformi ng cold gas into hot gas in the galaxies we dis- 
cusse d (e.g lMac Low k, Ferraral ll999; Du bois fc Tevssierl 
12008ft . 

5. CONCLUSIONS AND DISCUSSION 

We have investigated the recent star formation history 
of satellite galaxies by comparing semi-analytic models of 
galaxy formation with empirical data drawn from SDSS 
and GALEX. Based on the star formation r ate measure- 
ment s derived from multiband photometry ([Salim et al.1 
2007), we first divided galaxies into active and pas- 
sive types, and computed the fraction of passive galax- 
ies (/pass) as a function of gala xy stellar mass. As al- 
ready shown in previous studies ( Weinmann et al.l [20061 : 
iBaldrv et all 120061 : iKimm et all |2009j), satellite galax- 
ies in theoretical models cannot reproduce the observed 
level of recent star formation activity. The satellite 
over-quenching problem is generally attributed to the 
strong strangulation applied in most semi-analytic mod- 
els in which diffuse gas of the satellite system is instantly 
shock- heated from the system at the beginning of halo 
mergers (e.g. lWhite fc FrenklflOfll . 

In the hope of improving the situation, we implement 
gradual diffuse gas stripping and stellar mass loss into 
our semi-analytic code. For detailed tests, we divide the 
sample by morphology. We have also introduced an ap- 
proximation for the contribution to the ISM from stellar 
mass loss. Our results can be summarised as follows. 

• The over-quenched satellite galaxies in SAMs are 
mostly late-type. 

• The models with gradual diffuse gas stripping re- 
solve much of the satellite over- quenching prob- 
lem. However, they cannot account for the signif- 
icant fraction of actively star-forming small late- 
type galaxies observed in massive haloes. 

• Stellar mass loss is not a dominant source of cold 
gas in most satellite galaxies but enhances residual 
star formation. 

• Our new models incorporating both effects show a 
significantly-improved match to the observed data. 
However, they still suggest that the majority of 
cluster late types are passive. 

Recent SAMs including ours have shown steady progress 
in matching the observed data, but are missing key ingre- 
dients. For example, feedback from supernova explosions 
affects the recent star formation history in satellite galax- 
ies by regulating the remaining cold disc gas but is still 
poorly constrained (§4). 

One might wonder whether or not simply increasing 
the contribution of stellar mass loss to the cold ISM can 



reproduce the passive late-type fraction. We have found 
that such an appro ach still shows a lack of actively star- 
forming late types. iKavirai et al.l (|2007| ) reached a simi- 
lar conclusion when trying to explain early-t ype galaxies 
with residual star formation |Yi et al.ll2005f) . 

Recently, iTonnesen fc Bryanl (|2009l ) demonstrated 
that a weak ram pressure does not only remove cold gas 
in satellites but also may enhance star formation by com- 
pressing the gas component. Provided that the cooling 
of cold gas onto a galaxy disc is suppressed during the 
interaction between the host and satellite haloes and the 
satellites are affected by the weak ram pressure for a long 
time, this may increase / paS s by consuming a significant 
fraction of cold gas. Yet, it is still unclear how gas cool- 
ing takes place in massive satellite systems during halo 
merging. Since cooling also relies on the density, if ram 
pressure compresses both cold and diffuse gas, cooling 
could also occur more efficiently If this is the case, the 
efficient cold gas consumption might have little effect on 
/pass- Detailed numerical simulations are necessary to 
better understand the supply of cold gas during interac- 
tions. 

An accurate determination of initial orbits may have an 
impact on /p a s s . Based on numerical simulations at z ~ 
(iBensonl I200l iZentner et~aT1 [20051: iKhochfar fc Burkertl 
2006), we have assumed that satellite galaxies have ran- 
dom circularities at the time of halo mergers over the 
entire cosmic history. I f radial orbits are more commo n 
at higher redshifts (e.g. lWetzelll2010l : iDekel et al.ll2009f ). 
gas stripping due to ram pressure would be more effec- 
tive because satellites penetrate dense intra-cluster me- 
dia more frequently on radial orbits. Motivated by this 
idea, we performed a simple experiment: we assigned 
an eccentricity of I to the haloes that experience halo 
merger at z > 1. In later mergers, orbital parameters 
were randomly chosen in terms of circularity, as done 
in the fiducial model. Since massive galaxies are likely 
to orbit several times while smaller galaxies have passed 
the pericentre less, the influence of ram pressure on ra- 
dial orbits is more notable in massive satellite galaxies. 
As a result, the negative correlation of / paS s with m ga i 
shown by our new late-type galaxy models in massive 
haloes gets diminished, making the models match the 
data slightly better. However, the reliability of our exer- 
cise depends on the validity of the demarcation redshift 
and actual eccentricity distribution assumed. 

We also note that ignorance of the tidal disruption 
could have an i mpact on t he passive galaxy fraction . 
iTavlor fc Babul (f200l and IZentner fc Bullockl (|200l 
showed that subhaloes lose 30-40 per cent of their mass 
per pericentric passage. This implies subhaloes are dis- 
solved into host haloes after several o rbits unless they 
merge with centrals. On this basis, iSomerville et al.1 
(2008) implemented a prescription that satellites are dis- 
rupted when the subhaloes lose ~ 90% of their ini- 
tial mass. By doing so, they reproduced the luminos- 
ity and the radial d istri bution of Milky Way satellites 
(iMaccio et al.l I2010T ) . iHenriques. Bertone fc Thomas! 
(2008) also present good matches to galaxy colours using 
a simple assumption that satellites which are not associ- 
ated wit h subhaloes in the M illennium dark matter sim- 
ulation ([Springel et al.l 120051 ) are already disrupted by 
the tide in their envi r onme nt. Yet, as discussed in K09, 
the ISomerville et al.l (|2008l ) models still show a lack of 
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actively star forming small satellites, implying that ig- 
norance of tidal disruption is not the primary cause of 
red and dead small (late-type) satellites. Nevertheless, 
since it is small red and dead galaxies that are prefer- 
entiall y disrupted by tidal forces, the iSomerville et al.1 
(2008) models exhibit a slightly better agreement with 
empirical data than other semi-analytic models without 
tidal disruption (see K09 for details). In this regard, the 
passive galaxy fraction in the low-mass regime could be 
even smaller, resulting in a weaker dependence of the 
passive fraction on galaxy mass. 

Recent semi-analytic models often adopt AGN feed- 
back to prevent ho t gas cooling in massive galaxies (i.e. 
iCroton et al.ll2006l ). The effect has been particularly im- 
portant for "central" galaxies that are fuelled by ongoing 
cooling, in contrast to satellite galaxies that cannot re- 
tain their hot halo gas in most SAMs. However, the 
strong strangulation appears problematic in that even 
late types suffer from the lack of cold gas. For that rea- 
son, we demonstrate that external as well as internal gas 
supply are necessary to match the observed fraction of 
passive galaxies. Allowing the external supply has an 
interesting implication. Since the gas re tained in the 
hot halo could funnel into black holes (e.g ICroton et all 



2006), AGN feedback is likely triggered, and possibly 
suppresses star formation activity in satellite galaxies 
with supermassive black holes. On the other hand, the 
effect of AGN feedback may be negligible in late-type 
galaxies because their black holes are not massive enough 
to release a large amount of energy. Instead, the evolu- 
tion of late-type galaxies in clusters or groups may be 
more likely driven by the combination of stellar mass 
loss and environmental effects such as tidal and/or ram 
pressure stripping. We have indeed shown to first or- 
der that both mechanisms could reproduce observed lev- 
els of recent star formation, but it is still unresolved 
how late-type satellite galaxies show similar fractions of 
passive galaxies over various galaxy stellar masses for a 
given halo. Explaining these observations with realistic 
assumptions is an important challenge, especially since 
most of the galaxies in the universe are satellites. 
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